Thermalization of Interacting Fermions and Delocalization in Fock space 
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By means of exact diagonalization, we investigate the onset of 'eigenstate thermalization' and the 
crossover to ergodicity in a system of ID fermions with increasing interaction. We show that the 
fluctuations in the expectation values of the momentum distribution from eigenstate to eigenstate 
decrease with increasing coupling strength and system size. It turns out that these fluctuations 
are proportional to the inverse participation ratio of eigenstates represented in the Fock basis. We 
demonstrate that eigenstate thermalization should set in even for vanishingly small perturbations 
in the thermodynamic limit. 
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Introduction. - Statistical physics relies on the as- 
sumption that the system under investigation is in ther- 
mal equilibrium. However, what are the precise condi- 
tions for an isolated system to relax to thermal equi- 
librium? This question has a long history including 
the ground breaking numerical experiments initiated by 
Fermi, Pasta and Ulam p^ on an anharmonic chain of 
classical oscillators, where thermalization was not ob- 
served as expected |2. Nowadays, the investigation of 
thermalization in quantum many-body systems attracts 
a lot of theoretical attention, inspired by the new exper- 
imental possibilities in systems of cold atoms PHSJ- 

The trajectory of a classical ergodic system reaches 
all regions on the energy shell for sufficiently long times, 
establishing the microcanonical ensemble. As a conse- 
quence, suitable chosen subsystems obey the Boltzmann 
distribution. In the quantum case, switching on an inter- 
action in a many-body system will combine the unper- 
turbed eigenstates \i) of similar energies into new energy 
eigenstates: |a) ~ X]i"^f N)- If the expectation values 
Aa = {a\A\a) of observables in these new eigenstates 
approach their microcanonical values A-axicmiE) , as ob- 
tained by averaging over all unperturbed states in a small 
energy window around E, then the properties of ther- 
mal equilibrium are established in each many-body eigen- 
state. This is the essential idea behind the 'eigenstate 
thermalization hypothesis' (ETH) [5)17]. 

Recently, the ETH has been tested in numerical experi- 
ments [SUlOj. by means of exact diagonalization. For few- 
body observables like the momentum distribution, indeed 
it was demonstrated that Aa ~ ^micro(-E'Q) and that the 
fluctuations around A^icro decrease with increasing in- 
teraction strength and system size. 

In the present work, we address the important ques- 
tion of how fast thermal equilibrium is approached when 
increasing the system size. A direct, brute- force numer- 
ical approach would be prohibitive. Instead, we char- 
acterize the gradual delocalization of eigenstates in the 
many-body Fock space via the inverse participation ra- 
tio (IPR) Y.Zi(P?)^ (with P? = |cf n which turns out 
to be connected with the fluctuations of Aa- While a 
connection between the IPR and the fluctuations was 
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Figure 1. Probability distribution pf = l(Jja)P of a many- 
body eigenstate |a) in the non- interacting Fock basis (ener- 
gies ^i). a) For weak interaction Vjt = 0.1, the eigenstate 
is localized in Fock space, consisting of a few isolated peaks, 
b) At large V jt = 1.3, all Fock states with energies ^^ close 
to Ect contribute, c) pf averaged over a couple of nearby 
eigenstates in a range iEjt — 0.05, for V/t = 0.45, 1.45 (top, 
bottom). It can be approximated by a Lorentzian of width 
r (red line). The energy was chosen to correspond to infinite 
effective temperature (see main text). 



observed recently |1H [T2] , we are able to conjecture its 
functional form and its dependence on system size and in- 
teraction strength, based on earlier analytical results on 
Fock-space localization by P. Silvestrov. In particular, 
we have numerical evidence that the interaction strength 
needed for thermalization is below that needed for full 
quantum chaos. Moreover, we find that in the thermo- 
dynamic limit (TDL) thermalization (in the sense of the 
ETH) sets in for arbitrarily small interactions. This is in 
contrast to recent observations on relaxation in a classical 
ID system [13^. 

Here, we address these questions by means of exact 
numerical diagonalization for a system of spinless ID 
fermions on a lattice, where integrability is broken by 
an interaction of strength V. As the observable of in- 
terest, we consider the fermionic momentum distribution 
fk . Our main result is that for large enough V the fluctu- 
ations of /^ = {a\fk\c() are determined by the IPR which 
roughly can be considered as the inverse number of non- 
interacting Fock states \i) contributing to |a) (see e.g. 
|14|). The IPR itself keeps track of the transition from 
integrability to quantum chaos |11|[T5] and it was conjee- 



tured only recently that it might directly determine the 
deviations of steady state expectation values from the 
corresponding microcanonical value |16] . 
We observe three different regimes, depending on the in- 
teraction strength. An important scale is set by the 
mean level spacing A/ between Fock states that couple 
directly to a given initial Fock state. If the interaction 
is smaller than A/, then the eigenstates are 'localized' 
in Fock space |T71 [T5] and experience only a perturba- 
tive correction due to the interaction (see Fig.flk). For 
couplings beyond Af, the eigenstates delocalize and re- 
markably the IPR decreases exponentially with 1^ on a 
scale that depends on A/ . This scale essentially decreases 
polynomially in particle number and system size. There- 
fore, we expect the fluctuations of /^ to be suppressed 
to zero in the TDL even for vanishingly small interaction 
strength, establishing eigenstate thermalization of the 
considered observable. Increasing the interaction even 
further, eigenstates become chaotic (see Fig.fTJD) and the 
IPR as well as the fluctuations in /^ decrease as the in- 
verse many-body density of states as it was conjectured 
in iniEl- These results should apply rather generically 
to few-body observables diagonal in the eigenbasis of the 
unperturbed Hamiltonian. 

Model. - We consider n spinless ID fermions with pe- 
riodic boundary conditions on a lattice of TV sites and 
with a next-nearest neighbor interaction breaking the in- 
tegrability of the system. The Hamiltonian reads: 



JV 



i?o + V = -i^c|c,+i-|-H.c. 

+y^(ri,-l/2)(ri,+2-l/2). (1) 



The eigenstates \i) of Hq with ^^ = (*|^|*) are given by 
the Fock states of n fermions in momentum space. Due to 
the translational symmetry, the interaction does not mix 
Fock states with different total momentum K. Therefore, 
each momentum sector K with dimension Dk will be 
considered separately. We exclude the K — 0-sector as it 
possesses a trivial extra symmetry under reflection. In 
our numerical examples n — 7 and iV = 21. 

Fluctuations and IPR. - In the following, we discuss 
the expectation values f^ of the momentum occupation 

numbers fk = c\ck (where Cfc = l/\/^X]i=i ^ ''^^^^j). 
Being interested in the properties of typical eigenstates, 
we analyze the statistics of an ensemble of states \a) with 
similar eigenenergies Ea € Ie = [E ~SE,E + 6E], which 
will be called in the following 'eigenstate ensemble' (EE). 
The width of the energy window SE has to be chosen 
small enough to avoid artifacts resulting from systematic 
dependencies on E. Averages with respect to the EE are 
denoted by {...)e- For not too large interactions, one 
can easily show that {fk)E ~ fk,vaicm[E). However, the 
crucial statement of the ETH is that for each eigenstate 
itself f^ — )■ /fc, micro whcu going to the TDL, i.e., that the 
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Figure 2. Fluctuations of expected occupation number 
(aj/feja) between eigenstates decrease with increasing inter- 
action strength V, indicating eigenstate thermalization. Plot 
shows the variance Sfl^Q, for states with an effective temper- 
ature T/t = 1.2,1.7,00 (from top to bottom), with energy 
shells of width SE/t — 0.25. Solid lines show const x /j,^q(1 — 
fk^o) ^^ Var_B(pf ), with a slightly T-dependent constant. Fi- 
nally, we averaged the results over all total momentum sectors 
K. Inset: As in main figure, but with 5f^ averaged over all 
k. The red dots show Sfl averaged over all k for an inte- 
grable model with nearest-neighbor interactions (at T = oo 
and K/{2n/N) = 1). is'-averages are only performed to im- 
prove statistics. The same results are obtained for individual 
iC'-sectors. 



fluctuations of /^ from state to state vanish: 



Sfk 



{fk-fkY 



-> 0. 



(2) 



We introduced the EE-variance (5/| and fk = {fk)E- 
Representing f^ in the Fock basis /^ = J2i=iPifk 
(with fl = (i|/fe|i)) this statement becomes plausible. 
For strong interaction, typical eigenstates are spread out 
widely in Fock space (Fig.!!!)), i.e., they are composed 
of a large number of Fock states close in energy. Due 
to the law of large numbers, we thus expect the fluctua- 
tions to decay as the mean inverse number of Fock states 
contributing to \a), i.e., as the mean IPR 



Dk 
i=l 



(3) 



Before deriving the connection between J/| and x for- 
mally, we focus on the numerical results for the present 
model. Fig.p^ shows Sf"^ as a function of V evaluated 
w.r.t. eigenstates at various energies. The eigenener- 
gies can be re-expressed in terms of effective tempera- 
tures T, with Et = tiK{He-"/^)/tJ:K{e~"/^). The 
results are compared to the IPR, or more precisely to 



the sum over the variances Var£;(pf ) = {{pf)'^)E — {pf)^ 
(see discussion below), clearly demonstrating that indeed 
Sf^ oc X^i ^^''^^(pf ) Q'VG^ for small interactions. This is 
in stark contrast to the case of integrability conserving 
nearest-neighbor interaction (inset Fig. [21, where the sup- 
pression of (5/1 with V is much smaller than in the prior 
case. 

Formally, representing 6f^ in terms of pf, one finds 



<5/|^/fe(l-/,)^Vars(pn 



Dk 

E 



sfi^covEiPfpni) 



with Sfl'' — 
Covf 



iflfl 



k ~ fk ) E^nd the covariance matrix 
e{p?p'}) = {pfp'})E - {p?)E{pf)E- The first term 
in Eq. (W| contains the suppression of Sf^ with increas- 
ing number of Fock states contributing to a typical 
eigenstate. It is essentially determined by x [we note 
X ~ ^i^^^sipf) below the regime of full chaos (see 
below)]. We replaced ^.(/^ " 7')Vars(pf) ^ (/, - 

2 

fk)'l2i^^^E{pf), which is justified as Ya.rE{pf) is a 
smooth function of i. The pref actor fk{^ — fk) is nothing 
but the variance of the momentum occupation numbers 
for the non-interacting case. 

The off-diagonal contributions in Eq. (|4| are sensitive 
to residual correlations within eigenstates and are ex- 
pected to become small for strong perturbations. Sur- 
prisingly, for strong enough interactions, it approxi- 
mately reproduces the diagonal part of Eq. Q. Thus, 
even though 6f^ is still determined by the IPR, one ob- 
serves a deviation of the prefactor of 0(1). A very sim- 
ilar observation was made in [12] while investigating fi- 
nite fermionic systems with random two-body interac- 
tions and was traced back to the strong correlations be- 
tween matrix elements of two-body interaction matrices. 

To sum up, we find that the fluctuations in the ex- 
pectation value of fk from eigenstate to eigenstate are 
determined by the IPR x- Thus, in the following, it will 
be discussed how x decreases with increasing V and sys- 
tem size. Being a measure for the mean effective num- 
ber of Fock states forming an eigenstate, x indicates the 
'delocalization' crossover in Fock space and serves as an 
indicator for the transition from integrability to quantum 
chaos. 

Definitions - For the following discussion of the IPR, 
we need to set up a few technical definitions. We intro- 
duce the effective density p\{oj) of Fock states \j) cou- 
pling to a state \i) of energy ^^ G Ie (i.e., (J|V^|j) 7^ 0), 
where the energy difference between both states is ^i — 
^j = ut. Averaging over a couple of states \i) (indicated 
by (• • Ob = E^,4.e/E]"ME^,c.6/E ■ • ■]) one obtains the 
mean effective density of states pf{u!,E) — {p\{uj))^E. 
Furthermore, it will be convenient to introduce the in- 
teraction formfactor 
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Figure 3. a) The inverse participation ratio x decreases as 
many-body eigenstates get more delocalized for increasing in- 
teraction strength V. From top to bottom: Effective temper- 
atures T/t = 1.2, 1.4, 1.7, 2.6, CX3 (blue Unes; averaged over all 
momentum sectors). The most important feature is an ex- 
ponential decay at intermediate interactions (dashed line; see 
Eq. m) followed by a power-law tail in the quantum chaotic 
regime, where it can be well approximated by Eq. (|8| (shown 
only for T — 00 and K = 2-k/N: blue, dashed line). In the 
chaotic regime, the amplitudes cf are Gaussian distributed, 
leading to x ^ ^J2i{Pi)E (red lines), b) "Scaling plot": As 

before, but plotted vs. 2npfVV^ (and only for a single K). 
c) Similar plot, but for various system sizes, at T/t = 00. The 
red line displays Eq. ([6|, with C ~ 0.75, for comparison. 



This can be rewritten as i^ = npfV^, where V"^ denotes 
a mean matrix element squared. For w — > and small 
V, the form factor F reduces to Fermi's golden rule rate 
for a Fock state of energy E. In the following, only the 
mean matrix element and the effective density of states 
with respect to states close in energy, i.e., V'^{uj ~ 0,E) 
and Pf{uj ~ 0, E) will appear. For brevity these will now 
be denoted by V'^ and pf = A7^, respectively. 

Localized regime - As long as yV'^ ^ p7 , eigen- 
states can be obtained within standard perturbation the- 
ory (apart from a small set of eigenstates, which can be 
traced back to degenerate Fock states). A given Fock 
state gets perturbed by the set of directly coupling states 
and eigenstates consist of a small number of sharp peaks 
(Fig. Ilk), i.e., they are localized in Fock space. 

Delocalization - Increasing the coupling strength VV^ ~ 
p7 , one enters the regime of delocalized eigenstates |18| . 
Perturbation theory breaks down and the IPR starts to 
decrease rapidly (see Fig.lsk). In this regime, the fluctua- 
tions 5/| become directly determined by x- Surprisingly, 
one observes an exponential decay of x ^nd we found 
good numerical evidence that 



X(xexp{~CpfVV^}. 



(6) 



i=i;«7^i 



The numerical constant C is independent of temperature 
and system size. In Figs.|3]D,c, the IPR is shown as a 

function of the scaling variable 2TTpf\/V'^ for eigenstates 




Figure 4. a) Fock state decay rate F vs. V (at T = oo), ex- 
tracted from the imaginary part of the self-energy (red dots). 
For small V, Fermi's golden rule F = npfV^ holds, while F 
increases linearly in V for large V. For this plot, ImEi(tj) 
was averaged in both cj and energy ^i over the energy interval 
of width 5E — 0.25t centered around Et- Black dots show 
the results of a direct fit of {p")e- Here K = 2-k/N . b) Am- 
plitude distribution P(cf) for eigenstates at T/t = oo with 
5E/t = 0.25 and V/t = 1 demonstrating that for very large V 
one enters the chaotic regime. In this regime, the amplitudes 
cf are gaussian distributed as originally conjectured in iBj. In 
plane: energy dependent standard deviation ±[{(cf )'^)e]^ of 
P (blue lines). Out of plane: Cuts of P (black lines), which 
can be described by gaussians of variance {(cf )'^)b (red lines). 



at different energies E and for various N and n. Indeed, 
in good accordance to Eq. (|6| all curves collapse to the 
same scaling curve. An explanation of this exponential 
decay of x might be found in the two-particle nature of 
the interaction, following P. Silvestrov. In jlQj it was 
argued (in a random matrix setting) that for moderate 
interaction strength, typical eigenstates are composed of 
independent pairs of interacting fermions. Thus, eigen- 
states decompose into direct products of pairs of Fock 
states, resulting in an exponential decay of %, of the form 
given by Eq. (|6| . While this exponential decay (and ad- 
ditional corrections) have been confirmed numerically in 
a random quantum dot Hamiltonian [14] , here we find it 
in a translationally invariant many-body system without 
disorder. 

We now discuss the dependence on system size. The 
effective density of states p/(w) scales as N^. For ex- 
ample, at large T, we have Pf{uj)t ~ N^p^{l — p)'^r{uj), 
with the density p. For our particular model, r{uj) ex 
ln(i/u;) for w — >■ due to transitions of particle pairs 
around the infiection point of the — 2tcos(fc) dispersion, 
resulting in p7^ ex t/{N^lnN) (assuming a cutoff scale 
Lj/t ^ ^/N). Together with the scaling of the matrix 
elements \/v^{ll>,E) = v{uj,E)V/N, this would yield 
X oc exp(-CN^\nNp'^{l~ p)^V/t\, with C being in- 
dependent of iV. Thus, we expect the fluctuations to 
decrease drastically in the thermodynamic limit even in 
this intermediate regime, where eigenstates are not yet 
ergodic. 

Chaos - Only by increasing the interaction even fur- 
ther, one enters the regime of ergodic eigenstates. By 
'ergodic eigenstates', we understand states which in prin- 



ciple are composed of all Fock states close in energy (cf. 
Fig.IIlD). No Fock states are excluded a priori, e.g., due 
to the two-body nature of V or further symmetries from 
contributing to an ergodic eigenstate. The amplitudes cf 
become Gaussian distributed random variables [5117] as it 
is shown in Fig.l4|3 with a Lorentzian variance |HTJ [55] 



(P?) 



TiE~C.,E) 



npK{E)^^^_E-6{EM)?+T^' 



(7) 



where pK denotes the full many-body density of states 
for total momentum K , scaling as [N — 1)\/{N — n)\n\. 
This indicates the crossover to full quantum chaos. We 
checked that in this regime the nearest neighbor level 
spacing statistics agrees with the GOE-Wigner surmise, 
characteristic for GOE random matrix ensembles. Due 
to the Gaussian distribution for cf, one finds x — 
3Ea(pf)l; resulting in 



1 



X 



2t: T{Q,E)pK{Ey 



(8) 



which is in fairly good agreement with the numerical 
results in Fig. (Isk), demonstrating the suppression of 
(5/1 by the inverse many-body density of states as it 
was conjectured in [SI |7]. The mean spreading width 
F (Fig.Uk) can be extracted from the Fock state self- 
energy S by averaging — ImE(^i,a;) over ^i,a; € Ie- S 
is obtained from G(^i,w) = (i||u; -I- iO+ — H]^^\i) via 
G = [uj + iO+ — ^i — S]^^. Fig.fTp shows a comparison 
of {pf)E and a Lorentzian of width F extracted directly 
from — ImS. 

The important question remains, how the second 
crossover scale (governing the crossover from delocal- 
ized to ergodic eigenstates) depends on system size. We 
found some indication that it might depend on the in- 
tensive 'energy range' W of the coupling matrix V. Con- 
sider the dependence of F on V" in Fig.|4^. For small 
V, Fermi's golden rule applies and one finds r(0, E) ~ 
F{0,E) ex V'^/t. For large V, one observes a crossover 
r (X V^/t ^ T (X V indicating the entrance into the 
strong coupling regime, where F and the finite width (in 
w) of the formfactor F become comparable [2I4 . Com- 
paring Figs. [3^ and [4^, there might exist a close rela- 
tion between this crossover and the onset of ergodicity 
of eigenstates. This would imply that the interaction en- 
ergy scale pT for the onset of thermalization is paramet- 
rically smaller than the scale for the transition to chaos, 
determined by W. 

Conclusions. - By means of exact diagonalization we 
investigated the interaction induced onset of eigenstate 
thermalization in a system of ID fermions. We found 
that the fiuctuations of the expectation value of the mo- 
mentum occupation number from state to state are pro- 
portional to the inverse participation ratio of eigenstates. 
For small interactions the latter decays exponentially be- 
fore one enters the chaotic regime. The interaction scale 
for the onset of this decay is essentially set by the effec- 
tive mean level spacing between interacting Fock states. 



and this vanishes in the TDL. Thus, we corroborate the 
physical expectation that in the TDL at arbitrarily small 
interactions, eigenstate thermalization sets in. 
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